Dr. F.J. Rodenburg
June 2024 (rev. 1)
A new cancer treatment is compared to an existing alternative.
Example table of patients:
| patient | time | status | age | sex |
|---|---|---|---|---|
| 1 | 5 | dead | 65 | female |
| 2 | 11 | dead | 74 | male |
| 3 | 11 | dead | 81 | male |
| 4 | 11 | dead | 67 | male |
| 5 | 12 | dead | 74 | male |
| … | … | … | … | … |
Example table of patients:
| patient | time | status | age | sex |
|---|---|---|---|---|
| 1 | 5 | dead | 65 | female |
| 2 | 11 | dead | 74 | male |
| 3 | 11 | dead | 81 | male |
| 4 | 11 | dead | 67 | male |
| 5 | 12 | dead | 74 | male |
| … | … | … | … | … |
Example table of patients:
| patient | time | status | age | sex |
|---|---|---|---|---|
| 1 | 5 | dead | 65 | female |
| 2 | 11 | dead | 74 | male |
| 3 | 11 | dead | 81 | male |
| 4 | 11 | dead | 67 | male |
| 5 | 12 | dead | 74 | male |
| … | … | … | … | … |
Example table of patients:
| patient | time | status | age | sex |
|---|---|---|---|---|
| 1 | 5 | dead | 65 | female |
| 2 | 11 | dead | 74 | male |
| 3 | 11 | dead | 81 | male |
| 4 | 11 | dead | 67 | male |
| 5 | 12 | dead | 74 | male |
| … | … | … | … | … |
For each time point
Survival analysis is used to study time until event data.
library("survival")
require("ggsurvfit")
# Fit a Kaplan-Meijer curve
KM <- survfit2(Surv(time, status) ~ sex, data = lung)
# Plot the result
ggsurvfit(KM) + labs(x = "Days", y = "Probability of survival") +
add_risktable() + add_confidence_interval() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),
panel.background = element_blank(), axis.line = element_line(colour = "black"))
# Perform a log-rank test
survdiff(Surv(time, status) ~ sex, data = lung)
# Estimate survival at a given time (days)
summary(KM, times = 250)